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ABSTRACT 

In the fundamental quest of the rotation curve of the Milky Way, the tangent-point method has long been the simplest way to infer 
velocities for the inner low-latitude regions of the Galactic disk from observations of the gas component. In this article, we test the 
validity of the method on a realistic gas distribution and kinematics of the Milky Way, using a numerical simulation of the Galaxy. We 
show that the resulting velocity profile strongly deviates from the true rotation curve of the simulation because it overstimates it in the 
central regions and underestimates it around the bar corotation. In addition, its shape is strongly dependent on the orientation of the 
stellar bar with respect to the observer. The discrepancies are caused by the highly nonuniform nature of the azimuthal velocity field 
and by the systematic selection by the tangent-point method of high-velocity gas along the bar and spiral arms, or low-velocity gas 
in less dense regions. The velocity profile only agrees well with the rotation curve beyond corotation, far from massive asymmetric 
structures. Therefore the observed velocity profile of the Milky Way inferred by the tangent-point method is expected to be very close 
to the true Galactic rotation curve for 4.5 S K < 8 kpc. The gaseous curve is flat and consistent with rotation velocities of masers, red 
clump, and red giants stars measured with VLBI astrometry and infrared spectroscopy for R > 6 kpc. Another consequence is that 
the Galactic velocity profile for R < 4 - 4.5 kpc is very likely flawed by the nonuniform azimuthal velocities and does not represent 
the true Galactic rotation curve, but instead local motions. The real shape of the innermost rotation curve is probably shallower than 
previously thought. Using an incorrect rotation curve has a dramatic effect on the modeling of the mass distribution, in particular for 
the bulge component, whose derived enclosed mass within the central kpc and scale radius are, respectively, twice and half of the 
actual values. We therefore strongly argue against using terminal velocities or the velocity curve from the tangent-point method to 
model the mass distribution of the Milky Way. The quest to determine the innermost rotation curve of the Galaxy remains open. 
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1. Introduction 

The quest of the velocity field of the Milky Way has long been 
a very difficult task. The location of the solar system inside the 
disk makes it impossible to directly constrain the 3D position 
and velocity phase space for the whole Galaxy. Most surveys de¬ 
voted to the chemistry, kinematics, and dynamics of the Galactic 
disk(s) are based on optical spectroscopy, which is not appro¬ 
priate to probe the innermost low-latitude disk regions. It can 
only be surveyed at infrared, mm-, and cm-wavelengths where 
obscuration by dust is less problematic. This is one of the ob¬ 
jectives of the infrared Sloan Digital Sky Survey Ill’s Apache 
Point Observatory Ga lactic Evolutio n Experiment (APO GEE 
lEisenstein et alJ 1201 lib For instance, iBovv et all (1201 2l) have 
modeled the kinematics of 3365 red giants and red clump stars 
from the first year of APOGEE data and found a rotation curve 
consistent with a power-law model almost perfectly flat within 
4 $ R £ 14 kpc. Line-of-sight velocities of stars that yield indi¬ 
rect estimates of distances and rotation velocities are not efficient 
in probing the rotation curve in the innermost regions, however, 
which would be necessary to establish accurate mass distribu¬ 
tion models. The best technique currently in use to do this is 
probably astrometry of methanol and water masers in high-mass 
star formati on regions with Very Long Baseline Interferometry 


(VLBI, e.g., Brunthaler et al 

2009; Reid et al. 

2009; Asaki et al. 

2010) Shiozakietak|2011:1 

onma et al. 2012 

Rvgl et al. 2012 

Xu et all 2013; Zhang et al 

20141). Preliminary results from 


these experiments combined within the Bar and Spiral Structure 


Legacy (BeSSeL) Survey have for the first time directly re¬ 
vealed the spiral structure of the Galaxy and showed an appar¬ 
ently shallow inner rotation curve ( IReid et al.ll2014 ). Although 
very promising, it is too early for BeSSeL to entirely constrain 
the innermost rotation curve at a statistically significant level, 
however, since the current number of reliable velocities remains 
small, with fewer than ten sources inside R = 4 kpc. 

Instead, CO and Hi observations of the gaseous in¬ 
terstellar medium have long been used as the refer¬ 
ence to establish the inner rotation curve of the Milky 


1985; Fich et al. 

1989; McC 

ure-Griffiths & Dickev 2007; 

Levine et al. 2008!; Sofue et al. 

2009; Marasco & Fraternali 


20121) . Indirect estimates of rotation velocities and distances 


are obtained from the tangent-point (TP) method, which uses 
terminal line-of-sight velocitites (see Sectj2]i. Such veloci¬ 
ties and the rotation curve inferred from the TP method 
are t he kinematical basis of many Galactic mass models 
(e.g„ Merrifieldl 1 9921 Dehnen & Binnev 1998t IKalberla 20031 : 


iFamaev & Binnevl 20051 : ISofue et alJ 120091: Sofue 2013 ), 3D 
models of the inner neutral atomic and molecular gas disks of th e 
Galaxy (iNakanishi & So fue 2003, |2005 IKalberla et alJl2007h , 
or Galactic population synthesis models dRobin et al.ll2003l )~ 
However, this curve does not perfectly reflect the reality be¬ 
cause the gas distribution and kinematics are perturbed by the 
Galactic bar. The effects of the bar on gas terminal velocities 
and longit ude-velocity diagrams i n the central regions were firs t 
shown bv iLiszt & Burto n (1980) and iBlitz & S nergel ( 199 ill . 
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Fig. 1. Geometry of the system: Sun, tangent-point (T), Galactic 
Center (GC) assumed by the tangent-point method. The Sun is 
at radius Ro, T at radius R. The circular velocities are vo at Rq, 
and v at R , Z being the Galactic longitude. We choose here a 
view similar as in Renaud et al.l ( 2013 ) where the Northern and 
Southern Galactic Poles are respectively in front of and behind 
the reader. 




Fig. 2. Examples of spectral decomposition of the Hi profiles at 
Z = -9° and Z = 12°. The solid black curve shows the observed 
profile, which is decomposed using multiple Gaussian functions 
(in blue). The reconstructed profile is shown as a dashed red 
curve. The vertical dashed line shows the terminal velocity of 
the peak selected to derive the rotational velocity with Eq. [2] 


Numerical simulations of the Galactic interstellar medium have 
shown the dependency of such diagram^ on the position of the 
observer in the dis k dFuxl 19991: lEnglmaier & Gerhard 1999: 


iBissantz et alJ2003tlRodriguez-Fe rnandez & Combes 200§|). As 

a consequence, most of mass models restricted the analysis to 
|/| > 15 - 20° or R > 0.2Rq (Z is the longitude, R is the 
Galactocentric radius and Rq the Galactocentric radius of the 
Sun), to minimize the perturbing effects from the bar at the cen¬ 
ter. Although laudable, these arbitrary limits remain question¬ 
able because the bar and the spiral structure are expected to have 
a significant effect over more extended regions. 

In this paper we propose to test the capacity of the 
tangent-point method in predicting a reliable rotation curve us- 
ing_a hydrodynamical simulation of a Milky Way-like galaxy 
(iRenaud et alJl2013h . This simulation sets up a live model of the 
Galaxy in a self-consistent way, that is, without a fixed pattern 
speed. The resulting structure and velocity field of the stellar and 
gaseous components are far from being axisymmetric and circu¬ 
lar. In particular, noncircular velocities as high as 200 km s 1 are 
found along the major axis of the bar. We investigate the effect of 
the nonaxisymmetric perturbations on the shape and amplitude 
of the inferred curve and on mass models. We analyze the prop¬ 
erties of tangent points and constrain the radial range over which 
the tangent-point method provides a reliable estimate of the ac¬ 
tual rotation curve. The basics of the tangent-point method are 
described in Sect. [2] the confrontation to the gas simulation and 
a discussion of the limits of the method is provided in Sects. [3 
and[4] and the consequences for mass distribution modeling and 
for the Galaxy are detailed in Sect. [3 


2. The tangent-point method applied to the Galaxy 

For any given Galactocentric radius R smaller than the solar ra¬ 
dius, there are two line-of-sights at longitude +/ in the Galactic 
plane that are tangent to the circle of radius R (Fig. (TJ). The 
tangent-point method stipulates that the line-of-sight (1-o-s) ve¬ 


locity i' 7 -,ios is an extremum at the position of the tangent point 
T. The Galactocentric radius and circular velocity profile v(R) 
of such material at Galactic latitude b — 0° are expressed as a 
function of / and the terminal velocity v'r.ios by 


R — Rq sin(Z), 


( 1 ) 


v(R) = v 0 sin(Z) + v r ,i 0 s, 


( 2 ) 


where v ; o is the circular velocity at the Galactocentric radius Rq. 
It is straightforward to estimate the rotation velocity profile from 
Eq. |3 with mm- or cm- observations by selecting the Hi or CO 
spectral component whose 1-o-s velocity vr.ios is lowest for / < 0° 
and highest for / > 0°. 

We have applied the TP method to the Hi 
Leiden/Argenti ne/Bon n Survey datacube of the Milky Way 
(iKalberla et al.l 1200511 . This survey yields the most sensitive 
all-sky dataset of the neutral atomic gas of the Galaxy at angular 
and spectral samplings of 0.5° and 2 km s 1 . We used Rq = 8 kpc 
and vo = 213 km s 1 to facilitate comparison with the numerical 
simulation (Sect|3|. These are not the IAU values of 8.5 kpc 
and 220 km s _1 , but the choice of the fundamental parameters 
does not affect the result. Hi spectra were decomposed on a 
base of Gaussian functions, following the method used for the 
Andromeda galaxy dChemin et alJl2009h . In the example shown 
in Fig. [3 the spectra were fitted with nine Gaussian profiles. The 

selected velocity extrema of-205 km s~ ! at Z = -9° (R — 1.25 

kpc) and ~ 150 km s _1 at Z = 12° (R = 1.66 kpc) correspond to 
velocities of 238 km s -1 and 194 km s , respectively. 

By repeating this exercise for -90° < Z < 90°, we obtained 
the velocity curves shown in Fig. [3 which agree very well with 
many other curves using the same method from CO or Hi data 
fe.gJBurton & Gordonl 19781: IClemenslI 19851: ISofue et akl f2009: 


iMarasco & Fraternali l2012h . The steep profile in the center is 
clearly visible, as is the peak at ~ 260 km s' 1 at 500 pc, the 
smooth decrease to ~ 200 km s~\ the increase beyond R ~ 3.5 
kpc, and flat profiles at large radii. The axisymmetry between the 
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Fig. 3. Velocity profile of the Milky Way inferred with the 
tangent-point method applied to the Hi LAB datacube of 
dKalberla et al.l l2005h . assuming Ro = 8 kpc and vo = 213 
km s 1 . Blue upward triangles are for -90° < l < 0°, red down¬ 
ward triangles for 0° < l < 90°. The solid line is the composite 
curve from Sofue et al. ( 2009 ) normalized at our value of vq. 



Radius (kpc) 


Fig. 4. Rotation velocity profile of the simulated disk inferred 
by the tangent-point method for the axisymmetric disk (before 
the formation of the bar and spiral arms). Colored symbols are 
the same as in Fig. [3] The solid line is the average curve from 
the two halves. The dashed line is the true rotation curve of the 
simulated disk. 


two quadrants is quite remarkable for R > 3.5 kpc. The disk at 
l < 0° appears to rotate faster than at positive longitudes by ~ 17 
km s 1 on average for 1 < R < 3.5 kpc. Note in particular the 
bump at R ~ 1.5 kpc for l < 0° that is not measured for l > 0°. It 
has often been argued that differences in the central regions are 
signatures of perturbations, like the bar and/or the spiral struc¬ 
ture. This is a possible explanation, but we show in Sect. 13.21 
that a velocity difference is not the best criterion to detect sig¬ 
natures of perturbations. Notice finally that the composite curve 
of I Sofue et alJ (l2009t) is more consistent with our velocities for 
/ > 0° than for l < 0°. 

3. The tangent-point method to the test of 
numerical simulations 

3.1. Mock gas datacubes 

To quantify to which extent the TP method gives a reasonable 
rotation curve, we applied it to a simulation of the Milky Way, 
for which the rotation curve can be obtained directly from the 
velocity field. We used the simulation of iRenaud et al.l (1201 3i), 
which reproduces the main structures (bar, spirals) of the Galaxy 
both in terms of morphology and kinematics. We produced mock 
observations by computing longitude, latitude, and line-of-sight 
velocity datacubes for an arbitrary position of the Sun along 
the Ro = 8 kpc circle. Varying the position of the Sun allows 
us to monitor the effects of the bar orientation on the results 
from the TP method. In our fiducial case, the bar major axis 
yields an orientation of 23° with respect to the direction of the 
Galactic center. This angle setup matches the real Galaxy best 
(Renaud et al.ll2013h . The covered range of bar orientations is 0° 
(Sun aligned with the bar) to 90 ° (Sun - Galactic center axis 
perpendicular to the bar major axis). It allows us to test, for 
instance, whether the b ar viewing angle of ~45° inferred from 
some IR observations (Hammerslev et al. 2000; Beniamin et al.1 
1 20051) would yield results different from our fiducial case. The 
velocities are computed from the gaseous (atomic and molecu¬ 
lar) component of the model. The spatial and spectral samplings 
of these mock datacubes are 2° for (l,b) and 3 km s 1 for 1-o-s 
velocity. These samplings are smaller than for the Hi observa¬ 


tion, implying smoother velocity profiles, but this does not affect 
the result of this work. A datacube was also generated using an 
earlier epoch of the simulation. It corresponds to an axisymmet¬ 
ric disk, that is, a reference case before the formation of the bar 
and the spiral arms. 

3.2. Velocity profiles 

Figure |4] shows the resulting velocity profile for the reference 
datacube of the axisymmetric disk. The velocity increases within 
the inner 1.5 kpc and then remains constant with radius. The 
velocity curve averaged over positive and negative longitudes is 
exactly the one inferred from averaging the azimuthal velocity 
field of gas cells from the simulation at every radii. This implies 
that there are no systematics in our approach and that the TP 
method can yield a velocity profile that excellently agrees with 
the true rotation curve if the disk remains nearly axisymmetric. 

Using the snapshot that best matches the observations, that 
is, after the formation of the bar and the spiral arms, Fig.0shows 
the velocity profiles for a selection of several bar viewing an¬ 
gles. The shape and maximum amplitude of the velocity profile 
strongly depend on the bar viewing angle. The more aligned the 
reference point with the bar major axis, the higher the velocity 
peak at small radius. Moreover, the slope of the velocity profile 
in the central region decreases with increasing viewing angle. 
Configurations with angles between 0° and 23° all show a de¬ 
crease of velocity from R = 0.5 - 1 kpc to R ~ 4 kpc, while 
the velocity increases up to R = 3-3.5 kpc for larger viewing 
angles. The velocity peak is thus offset for large viewing angles. 
With the exception of the profile for an angle of 90° that continu¬ 
ously overestimates the true values beyond 3.5 kpc, most config¬ 
urations yield velocities that agree fairly well with the correct ro¬ 
tation curve for R > 4.5 kpc. In addition, the velocity is constant 
beyond 5.5 kpc. A common feature to most velocity profiles is 
the inflection of the curve past the bar corotation ( R c = 3.6 kpc). 

The velocity difference between the two quadrants strongly 
varies with the bar viewing angle. For instance, the two quad¬ 
rants yield very similar velocity curves at every radius for con¬ 
figurations of 15° or 45°, while larger differences are found at 
0° and 67°. The asymmetry is also stronger for R < 4.5 kpc 
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Fig. 5. Rotation velocity profile of the simulated disk for several viewing angles of the bar. The vertical dashed line is the radius of 
the bar corotation, as computed in lRenaud et al.l ( 2013 ). 


at 23° and 67°. The bar and spiral arms cannot be detected 
when viewed at 15° and 45°, with the assumption that differ¬ 
ences between the two quadrants trace nonaxisymmetric struc¬ 
tures. Similarly, it is difficult to reconcile the asymmetry at 0° 
and 67° at large radii where the perturbations are weaker with 
the axisymmetry at similar radii for all other configurations. 
Differences between the two quadrants thus cannot be the sole 
criterion to attest or reject the presence of nonaxisymmetric 
structures like the bar, or assess its strength. 

The comparison between the rotation curve of the mock 
Galactic disk and the velocity profiles from the TP method all 
given at (Ro, vo)=(8 kpc, 213 km s _1 ) is very instructive (Fig. ©. 
The tangent-point method rarely gives a velocity profile that 
agrees with the rotation curve for R < 4.5 kpc. The profiles for 
bar viewing angles 45° overestimate the rotation curve inside 
R = 2.5-3 kpc and then underestimate it past this radius. Almost 
all profiles are consistent with the true rotation beyond R = 4.5 
kpc. The extreme viewing angle of 90° is the only model that 
uniformly overstimates the rotation curve at these radii, while 
it understimates it in the inner regions. To illustrate the signifi¬ 
cant discrepancy at small radius, the difference of velocities for 
the fiducial bar orientation of 23° is 100% of the true velocity at 
R = 0.5 kpc, ~50% at R — 1 kpc, and 18% at R = 3.7 kpc, close 
to the bar corotation. 

Figure [6] also shows the comparison with the observational 
Hi data. The agreement is very good at R > 4 - 4.5 kpc for 
most of bar orientations. Bar viewing angles of 23° and 30° are 
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Fig. 6. Comparison of the true rotation curve of the simulated 
disk (filled symbols) with the velocity profiles inferred by the 
tangent-point method (colored lines). The viewing angle of the 
bar is indicated for each velocity profile. The open circles repre¬ 
sent the velocities from the Hi observations of the Milky Way. 


the models that best reproduce the central slope and the ampli¬ 
tude of the peak at R = 0.5 kpc. Smaller angles overestimate 
the peak, while larger angles underestimate it. Angles of 67° and 
90° are the least compatible models with the observed velocities. 
On average, the shape and amplitude of velocities in the fourth 
quadrant (/ < 0°, upper part of the observational points) are more 
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Fig. 7. Left: Gas density map showing the locations of terminal 1-o-s velocity points (white circles). The black circle indicates the 
positions of the tangent points. The reference position is (x, y) = (-8,0) kpc. Right: Comparison of the distance (top) and velocity 
(bottom) of the tangent points with the true distance and velocity of the terminal velocity points. A d = d-y - d where d-y — (Ry-R 2 ) i/2 
is the distance to tangent points and d the true distance to terminal velocity points. Av = v(R) - vj, , where v(R) is from Eq.[2]and 
is the true azimuthal velocity of terminal velocity points. 


consistent with the models than in the first quadrant (lower part 
of the points). 

4. Limits of the tangent-point method 

4. 1. Properties of tangent points and terminal velocity points 

Although the tangent-point method is appropriate for the early 
stage snapshot (Fig. [4}, its validity is clearly questionable for the 
asymmetric disk. To address this problem, we analyzed the char¬ 
acteristics of the tangent points and terminal velocity points. We 
based the analysis only on the fiducial case (bar at 23°) because 
it represents the real Galaxy best (Renaud et al.l2013t) . 

We first address the question whether tangent points are ac¬ 
tual terminal velocity points. This goal can be achieved by com¬ 
paring the tangent points with the positions and azimuthal ve¬ 
locities of the terminal velocity points. We find the gas at the 
origin of the terminal velocities by searching the clouds whose 
1-o-s velocity is within +5 km s from the terminal velocity. 
This velocity range corresponds to a window of 5cr, where cr 
is the spectral sampling of the mock datacube. Since more than 
one cloud may contribute to shape the emission profile around 
the terminal velocity, we expect to find a few lines of sight that 
harbor several clouds within the search window. In this case, we 
kept the cloud with the brightest contribution to the emission 
profile because it represents that population of clouds best. Note 
that a wider search window only marginally changes the identi¬ 
fied clouds. 

Figure [ 7 ] shows the corresponding locations of the terminal 
velocity points. These points are almost exclusively clouds lying 
along the spiral arms and the bar. A few positions are not asso¬ 
ciated with a specific cloud, but correspond to less dense gas in 
inter-arm regions (e.g., x,y ~ -1.3,2.5), or along a more dif¬ 
fuse spiral structure (e.g., x < -6, y > 0). However, it is striking 

1 We roughly estimate the brightness for the simulation as the gas 

density divided by the square of the distance to the Sun. 


that terminal velocity points that coincide perfectly with the tan¬ 
gent points whose locations are delineated by a circle are rare. 
In particular, the pitch angle of the spiral arms passing at the 
Sun position causes the tangent-point method to systematically 
overestimate the true distance at high longitudes. 

Figure [7] shows the errors on the true distance to the termi¬ 
nal velocity points ( d ) when assuming that the distance is that of 
the tangent points d-y = ( R 1 2 - R 2 f ' 2 (see Fig.Q}. On average, 
the difference in distance is 13% for -55° < l < 0° and 8% for 
0 ° < l < 55°. That ratio increases at |/| > 55°, where it is ~ 86% 
on average. Therefore, the closer to the Sun the tangent point, the 
less accurate the distance dy. Figure|7]also shows the difference 
between the velocity of the tangent point as given in Eq. [2] and 
the local azimuthal velocity at the position of the terminal ve¬ 
locity point. The tangent-point velocity underestimates (oversti- 
mates) the azimuthal velocity at negative (positive) longitudes. 
The difference increases toward low longitudes. On average, a 
difference of 7% is measured for -50° < l < -10° relative to the 
local azimuthal velocity; this is 4% for 10° < / < 50°. The error 
remains relatively constant within this longitude range. It is neg¬ 
ligible at high (absolute) longitudes and can be larger than 15% 
in the inner |/| < 10°. We conclude that the tangent points and 
terminal velocity points are thus quite similar for intermediate 
longitudes, but are certainly not the same at the lowest and high¬ 
est longitudes, or at particular outlier directions. The effect of 
the significant difference of position relative to the true position 
at large / is negligible for velocities, however, since the rotation 
velocity does not vary significantly at large radii. 

Similarly, we analyzed the differences between the tangent- 
point velocity and the local azimuthal velocity at the position 
of the tangent-points. Here we again find a strong dependency 
on longitude, but the scatter is about five times larger than what 
is shown in Fig. [7] for terminal velocity points. In other words, 
a tangent-point velocity is surprisingly less consistent with its 
local azimuthal velocity than with the azimuthal velocity at the 
position of the terminal velocity point. 
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Interestingly, we note that velocity profiles based either on 
the local azimuthal velocity (and position) at the terminal veloc¬ 
ity points or on the local azimuthal velocity at the tangent-points 
differ significantly from the rotation curve for R < 5 kpc. Their 
shape remains similar to v(R) of Figs.0and 0 again with veloc¬ 
ities overestimated in the inner regions, and is underestimated 
around the bar corotation. 

The distance errors have previously been estimated by 
iGomezI d2006l) . but from more idealized gas flow simulations 
than the one we used, for a d isk perturbed by a m — 2 spi¬ 
ral within the fixed potential of [Dehnen & Binnevl d 1998b . They 
have found mean errors lower than 0.5 kpc, and larger errors in 
the spiral arms. Our result is thus quite consistent with their es¬ 
timate, on average. 

4.2. Origin of the discrepancy with the rotation curve 

We have shown that Eqs. Q] and [2] are rough - and sometimes 
incorrect - approximations of the local positions and azimuthal 
velocities of terminal velocity points and are even poorer esti¬ 
mates of the local azimuthal velocities at the coordinates of tan¬ 
gent points. However, the position and velocity errors are not 
issues as such. Indeed, it remains difficult to recover the shape 
of the rotation curve of the mock disk, irrespective of the nature 
of the points or velocities chosen to calculate the velocity profile. 
At the same time, we have shown that the bar orientation with 
respect to an observer strongly affects the shape of the inferred 
velocity profile. Ignoring the effect of the bar therefore renders 
the TP method useless in regions where the bar dynamics dom¬ 
inates the motions, and one may thus accordingly question its 
validity in estimating the true rotation curve of a barred spiral 
galaxy like the Milky Way. 

The discrepancy with the rotation curve is not directly linked 
to the tangent-point method itself. It is a combination of several 
elements: the asymmetric nature of the disk, the nonazimuthal 
dependency of any rotation curve, and the location of the Sun. 
It is indeed the structure of the velocity field itself that causes 
the peculiar shape of the velocity profile: the tangential veloci¬ 
ties of gas are highly nonuniform. Perfect axisymmetry and uni¬ 
form tangential motions do not exist. This is evidenced by the 
strong dependency of tangential velocities on azimuth (Fig. [8j 
at R — 1 kpc, e.g.). The axisymmetric value is ~ 175 km s' 1 , 
whereas the bulk of gas rotates with velocities ranging from 100 
km s' 1 to 290 km s -1 , with a few points down to 50 km s' 1 . The 
rotation curve smears out such nonuniform tangential motions 
because a rotation curve only depends on R. With the tangent- 
point method, the velocity field is unavoidably restricted to two 
velocities at each radius, arising from the two quadrants. We 
showed that this is sufficient to derive the rotation curve for the 
early stage in the simulation since the tangential velocity is uni¬ 
form at that point, the disk being almost perfectly axisymmetric. 
However, this is not the case anymore for the barred and spiral 
configuration. The location of the Sun reference point (thus the 
viewing angle of the bar) causes the two velocities to coincide 
with high-velocity clouds along the bar and in the spiral arms, 
or low-velocity gas in lower density regions. They can therefore 
only be representative of local nonuniform tangential motions, 
not of the global rotation curve. Figure [8] illustrates this argu¬ 
ment at R = 1 kpc, where the velocities of the two terminal- 
velocity clouds coincide with the upper envelop of the curve. 
Larger (smaller) bar viewing angles than 23° would shift the tri¬ 
angles to the left (right) along the curve. At viewing angles of 
~ 67°, the velocities would be almost coincident with the hori- 
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Fig. 8. Top: Azimuthal profile of tangential velocity at R = 1 
kpc. Colored symbols are the velocities inferred by the tangent- 
point method at a similar radius for the two quadrants. The hor¬ 
izontal dashed line is the azimuthal average of the velocity pro¬ 
file, hence the velocity of the rotation curve at R = 1 kpc. The 
zero azimuth corresponds to the direction of the Sun. Bottom: 
Profile of radial velocity v«, locally given at the positions of the 
tangent points. Each value is not an azimuthal average from the 
radial velocity field at each radius. Symbols are the same as in 

Fig- El 


zontal line (azimuths of ~ 50° and 230°), and the rotation curve 
would be correctly recovered (see also Fig. 01. 

Confirmation of the local motions affected by the bar at low 
radius is also seen in the radial velocity vr profile, where vr are 
given here at the positions of the tangent-points (Fig.0 bottom 
panel). Significant noncircular motions down, for example, to 

-100 km s' 1 are detected exactly where the velocity profile of 

Fig.0differs the most from the rotation curve, as well as at larger 
radius. This is not fortuitous, as the asymmetries affect both the 
circular and noncircular motions. 


5. Discussion 

5.1. Implication for the mass distribution 

A straightforward implication of using the incorrect rotation 
curve based on the tangent-point measurements is the incorrect 
modeling of the mass distribution at the center of the galaxy. 
Many mass models of the Milky Way invoke a dominant con¬ 
tribution of the stellar bulge in the central regions to explain 
the steepness and the peak of the profile seen i n Fig. E] (e.g., 
iMerrifieldl 1 19921 ISofue et ail 120091 : ISofuel 1201 3l) . We proceed 
similarly with the simulation using the incorrect rotation curve 
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Fig. 9. Top: Mass distribution models with the true and incorrect 
rotation curve of the simulated galaxy. Open symbols and the 
green dashed line denote the true rotation curve and the bulge 
contribution, filled symbols and a solid line represent the incor¬ 
rect rotation curve and bulge contribution. For clarity, the re¬ 
maining contribution from the gaseous and stellar disks plus dark 
matter is not shown. Bottom: Error on the integrated mass of the 
bulge. AM = Mi nc - M true , where M; nc is the mass deduced from 
the incorrect bulge, and M true the true bulge mass. 


from the fiducial case and fit a bulge that best reproduces the 
velocities for R < 0.8 kpc (Fig. [9]). 

The differences between the contributions from true bulge of 
the simulation and the incorrect fitted bulge are shown in Fig. [9] 
The true bulge contribution is not a fit to the true rotation curve, 
but comes directly from the identification of stellar bulge parti¬ 
cles in the simulation. 

The true and false bulges are modeled by a spherical density 
profile given by 

p(r) = P —2 exp (-2 n ((r/r_ 2 ) 1/n - 1)) . (3) 

It is equivalent to an Einasto model (iNavarro et al Jl2004h . where 
2 is the scale radius at which the density profile has a loga¬ 
rithmic slope of -2, p -2 is the scale density at that radius, and 
n a dimensionless index describing the shape of p(r). At fixed 
scale density and radius, the smaller the index, the shallower the 
density profile in the center ( iChemin et alJhOlTh . 

The parameters we found for the incorrect bulge are p_ 2 = 
10.15 10 _1 M e pW 3 , r_ 2 = 0.42 kpc, and n = 0.95. The veloc¬ 
ity peak of the incorrect bulge is 260 km s" 1 at R — 0.7 kpc, 
thus 54% higher than the true bulge peak, and at a radius twice 
smaller. Note that the velocity shape and amplitude of that in¬ 
correct bulge are comparable with those seen in mass models of 
the Milky Way ( ISofue et~akll2009t) . The Einasto parameters of 
the true bulge are p 2 = 9.46 10 1 M e pc~ 3 , r_ 2 = 0.87 kpc, and 


n = 1.14. The incorrect bulge is thus about twice smaller and 
slightly denser than the true bulge. As expected, its mass dis¬ 
tribution is therefore more concentrated than the true one for a 
similar enclosed masses (~ 1.6 x 10 ll) M e for the true bulge and 
~ 1.7 X 10 10 M 0 for the incorrect bulge inside R = 8 kpc). The 
effect of this higher mass concentration on the enclosed bulge 
mass is shown in the bottom panel of Fig. [9] As rule of thumb, 
the considerable error on the enclosed mass is ~330% of the true 
mass inside R = 0.5 kpc, or ~ 140% inside R = 1 kpc, the in¬ 
correct enclosed bulge mass being 4.3 and 2.4 times higher than 
the true one inside these radii. Note that the stellar disk plus dark 
matter halo contribution remains similar in both mass models. 

5.2. Implications for the Milky Way 

In view of the results, the tangent-point method and terminal ve¬ 
locities very probably yield an incorrect rotation curve of the 
Milky Way inside R — 4 - 4.5 kpc. The steep central gradient, 
the peak at R = 0.5 kpc, and the smooth decrease to R = 3.5 
kpc reflect the highly nonuniform nature of tangential velocities, 
where gas locally orbits faster or slower inside the Galactic bar, 
the spiral arms, and the interarm regions than the average at the 
same radii. These local features cannot be representative of the 
true inner rotation curve of the Galaxy. Beyond R = 4 - 4.5 kpc 
(R/Ro > 0.5 - 0.56), it is reasonable to assume that the Galactic 
rotation curve is consistent with the velocity profile obtained by 
the tangent-point method, unless stronger asymmetries than ex¬ 
pected by the simulation exist at these radii in the Galactic disk. 
As a consequence, mass distribution models of the Milky Way 
based on terminal velocities in the inner regions are unavoidably 
flawed. We estimate that the Galactic bulge might be less con¬ 
centrated than currently thought, with a characteristic size of ~ 1 
kpc, thus about two times larger than inferred by mass models 
based on terminal velocities. 

Furthermore, the tangent-point method gives rise to an inner 
steep velocity profile and a peak under particular circumstances, 
when the bar is viewed with angles < 45°. The central velocity 
peak is always overestimated at viewing angles ^ 20°. Mock 
velocity profiles with bar orientations of 23° and 30° are the 
models that agree best with the observational velocities. Galactic 
bar viewing angles of ~ 23 - 30° are thus favored by our anal¬ 
ysis. This agrees with many observational or numerical studies 
(iRodriguez-Fernandez & Combesl2008UWegg & Gerhardl2013b 

Rena ud et al .112013l). but disagrees w i th the obs ervational value 
of ~ 44° found bv lHammerslev et al.l (|2000) and lBeniamin et afi 
(2005). In addition, the inflection of the observed velocity curve 
atR = 3.5 ±0.3 kpc probably marks the position of the corotation 
of the Galactic bar. 

We can deduce a hypothetical, but more realistic, Galactic 
rotation curve by naively applying a correction factor of ( 7 ^) . 
to the observed profile, which we deduce from the simulation 
(see Sect. l3.2l and Fig. [9}. This first-order approach assumes that 
the effects inducing errors in the estimate of the rotation curve 
from the TP method in the simulations are the same as in the 
real Galaxy. We have no direct way to verify the validity of this 
assumption and thus, the following hypothetical curve should be 
considered with caution. It gives for instance velocities that are 

2.2, 1.5 and 1.1 times lower than expected by the TP method at 
R = 0.5,1, and 2 kpc, respectively, and about 1.2 times higher at 
R = 3.5 kpc (for vo = 213 km s _I ). 

Figure [TO] compares individual rotation velocities from 
BeSSeL, the astrometric survey of methanol and water m asers 
in high-mass star-forming regions (Reid et al.ll2009i. 12014 ). and 


7 

























Laurent Chemin et al.: The incorrect rotation curve of the Milky Way 



Radius (kpc) 


Fig. 10. Hypothetical gaseous rotation curve of the Milky Way 
(solid line) after correcting for the initial velocity profile of 
Fig. [3] The cyan, pink, and gray shaded areas are the 1, 2, and 3<x 
confidence levels . Ope n circles are rotation velocities of masers 
from Reid et al] (l2014l> . The dashed line is the power- law rota¬ 
tion curve model of the Galaxy from lBovv et alJ (1201 2l) . 


the R > 4 kpc stellar rotation curve from APOGEE (Bow et al .1 
[2012), with our hypothetical gaseous rotation curve. Both stellar 
and gaseous velocities were rescaled here using vo = 218 km s _1 , 
as given in iBovv et al.l (i2012l) . The velocity uncertainty for our 
gaseous rotation curve takes into account the formal error from 
the Gaussian fits of terminal velocities (< 1 km s _1 ), and the sys¬ 
tematic errors from (i) the mean velocity difference between pro¬ 
files from simulated datacubes with bar orientations of 23° and 
30° (5 km s _1 ), (ii) the velocity difference Av = v(R) —v^ implied 
by the use of TP method velocities instead of true azimuthal val¬ 
ues and (iii) the velocity difference from both quadrants. These 
two latter systematic errors depend on the longitude, as seen in 
Figs. [3] and [7] (right panel). First, it is interesting to note that all 
rotation curves are consistent for R > 4.5 kpc, which clearly 
shows that the Milky Way rotation curve is consistent with a 
plateau at R > 6 kpc. Secondly, our hypothetical rotation curve 
tends to overestimate the BeSSeL velocities at 3.5 < R < 4.5 
kpc, though remaining in agreement within the 3cr confidence 
level. This is probably due to systematic effects inherent to the 
construction of our hypothetical curve, that is, to the TP method 
and the simulation themselves. Note, however, that depending 
on the location of the masers down- or up-stream of the peak 
densities inside the spiral arms, velocities from BeSSeL may not 
coincide perfectly with the mean rotation velocities either. They 
might trace faster or slower local motions, exactly like those seen 
for individual gas clouds selected by the TP method (Fig [8j see 
also Fig. 4 in lKawata et al.ll2014f) . It appears essential here that 
the inner Galactic velocity field has to be covered as uniformly 
as possible in radius and azimuthal angle to obtain the most ac¬ 
curate inner rotation curve, similarly as is routinely done for 
nearby galaxies with 3D spectroscopy or radio interferometry 


nearby galaxies with jD spectroscopy or radio interferometry 
( Emsellern^alJ 2004; Chemin et al.ll2006t Ide Blok et alJ l2008: 


iGarcfa-Lorenzo et alJl2015h . This would require observations of 

other tracers than masers in more diffuse regions outside the bar 
and the arms. This is not yet possible for VLBI, however. 

A final consequence of the analysis is that rotation curves of 
edge-on gala xies based on envelope-tracing of Hi terminal ve- 
locities (e.g JSwaters et al] fl997: Kre gel & van der Kruid [20041 
might be systematically biased toward too high or low velocities 
as well. Unless no bar or strong perturbation is shown to exist 


in their central regions, caution should be taken when model¬ 
ing these velocity curves because a steep (shallow, respectively) 
velocity rise could mimic a stellar component that is too (less) 
concentr ated, or a dark matter halo that is too cuspy (shallow; 
see also iRhee et aTll2004t I Valenzuela et al.l 120071 : iDicaire et al.l 
2008, for less inclined disk cases). 
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